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Abstract In this article we find all possible matrix means which are points of 
geodesies of affinely connected manifolds. We characterize certain properties of 
these manifolds, decide whether they are metrizable or not. We show that the 
matrix means that occur in this form are exactly the matrix power means. We 
, show that the several variable formulas of power means are the unique solutions 

of certain Karcher equations corresponding to affinely connected geometric struc- 
tures. At the same time we construct weighted matrix means corresponding to a 
matrix mean and we show that the representing operator monotone functions of 
the weighted means form a one parameter family of functions that has a semigroup 
structure in the sense of Loewner. These one parameter families are induced by 
operator monotone functions that behave similarly as logarithm maps of affinely 
connected manifolds. We show that the holomorphicity of the members of these 
one-parameter families is controlled by the distribution of the branch points of the 
corresponding logarithm map. 
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1 Introduction 

The open convex cone of positive definite n-by-n matrices, denoted by P(n, C), 
as a smooth manifold, carries some important Riemannian structures which are, 
in fact, symmetric spaces. For example consider the classical Riemannian sym- 
metric space P(n, C) = GL(n, C)/U(n, C) which plays an important role in several 
constructions [9J. This symmetric space is nonpositively curved, hence a unique 
minimizing geodesic between any two points exists. The midpoint operation on 
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this space, which is defined as taking middle points of the geodesic connecting two 
points, is the geometric mean of two positive definite matrices [6]. 

The other two symmetric spaces which are well known are Euclidean spaces. 
One of them is just the vector space of squared complex matrices. The subset 
P(n, C) therefore inherits its Euclidean structure and the midpoint operation is 
the arithmetic mean. The third case is the harmonic mean, the corresponding 
space is isometric to the previous one. 

These afRnely connected manifold structures are important and widely used 
in several cases. For instance the extension of the two variable geometric mean 
to several variables are heavily based on the Riemannian structure of the corre- 
sponding space [1,28,25]. The multivariable geometric mean or Karcher mean is 
defined as the center of mass 

n 

A{A 1 ,...,A n )= argmin V<5 2 (X,^). (1) 

on the Riemannian manifold P(n, C) endowed with the trace metric 

d(A,B) = y/Tr\og{A~ 1 B). 

This mean and the geometric structure corresponding to it is also used in practical 
applications when one considers averaging of symmetric tensors [2"l lllll2Tl |3l|3] . The 
Karcher mean A(A±, . . . , A n ) is also the unique positive definite matrix where the 
gradient of the function in the minimization problem (Q]) vanishes 

n 

J2^og(X- 1 A l ) = 0. 

i=i 

This equation is called the Karcher equation. Recently Lim and Palfia found 
a one parameter family of multivariable matrix means called the matrix power 
means which are defined as the unique positive definite solution of the equation 

1 - 

i=l 

where G t {A,B) = A 1 / 2 (AT 1 / 2 BA~ 1 l 2 ^j A 1 / 2 is the weighted geometric mean of 

A, B € P(n, C). The attractive property of this family is that as the defining 
parameter t — > 0, the matrix power means converge to the Karcher mean. In this 
paper in Section 7 we will show that the matrix power means are actually unique 
solutions of Karcher equations 

n 
i=l 

corresponding to affine connections given on the manifold P(n, C) . In other words 
the two-variable matrix power means give back the geodesic lines corresponding 
to these affine connections. Moreover in Section 8 we find all possible two-variable 
matrix means which occur such a way, in other words as affine matrix means, 
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which question was raised in [32]. It turns out that these means are exactly the 
matrix power means. We prove also that the corresponding affine connections are 

X7 Xp Y p = DY\p][X p ) - | (Xpp-Wp + Ypp^Xp) , 

where < re < 2 and the tangent space is H(n,C), the space of n-by-n hermitian 
matrices, at every point p € P(n, C). These connections appear earlier when we 
construct them as prototypes of invariant affine connections in Section 6. In Section 
9 among other results we show that these affine connections are non-metric in 
general, i.e. there exist no other Riemannian structures as in the case of the Karcher 
mean ([I]). In order to achieve this we investigate the holonomy groups and other 
properties of these affine connections. 

During the above process we give some geometric construction which can be 
used over a general affinely connected space to reconstruct the logarithm (hence 
also the exponential) map of the corresponding affine connection from the midpoint 
operation m(p,q) = exp p (l/21og p (g)) on the manifold as 

! / s r ■m(p,q)° n -p 
lo g p W = llm — 1 

2™ 

where we use the notation m(p,q) on = m (p, m(p, q) ^ 1 ^ 1 ^ and \og p (q) is the log- 
arithm map. This is dicussed in Section 3. We apply an analogue of such a process 
to arbitrary two- variable matrix means in Section 4 and we obtain a corresponding 
logarithm map 

log A (B) = A 1 ' 2 ^ U-V 2 BA-V*\ A 1 ' 2 

to the matrix mean, where \ogj(x) is an operator monotone function. We show 
that \og A (B), hence \og I (x) directly induce a one parameter family of operator 
monotone functions which represent matrix means that can be thought of the 
weighted counterparts of the original mean. Then in Section 5 we show that tak- 
ing directly such operator monotone functions \ogj(x) that can be prototype of 
logarithm maps, we obtain similar one parameter families of matrix means. We 
show that in both cases these families are Loewner semigroups of Pick functions 
which itself has a classical and rich theory [12, 13, 17, 18, 19, 20 . We also show that 
the further extendability to greater parameter values of the one parameter family 
depends on the distribution of the branch points of the corresponding logarithm 
map. We relate this extendability property to functional equations over the upper 
complex half-plane of the form 

lo gl (f t (z)) = t\o gl (z), (3) 

where ft(z) is the representing operator monotone function of the matrix mean 
and log^z) is the corresponding logarithm map. We show that if log/(z) has no 
branch points in the uper half-plane, then the functional equation, hence the the 
one parameter family ft(z) is a Pick function, i.e. an operator monotone function 
for all t e [0,1]. 
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2 Matrix means and some constructions 

Let us recall the family of matrix means [52] : 

Definition 1 A two-variable function M: P(n,C) x P(n, C) h-> P(n, C) is called a 
matrix mean if 

(i) M(I, I) = I where I denotes the identity, 

(ii) if A < A' and B < B' , then M(A, B) < M(A',B'), 

(iii) CM {A, B)C < M(CAC, CBC) for all hermitian C, 

(iv) if A n I A and B n I B then M(A n , B n ) I M(A, B). 

In property (ii), (iii), (iv) the partial order being used is the positive definite order, 
i.e. A < B if and only if B — A is positive semidefinite. An important consequence 
of these properties is |24| that every matrix mean can be uniquely represented by 
a normalized, operator monotone function f(t) in the following form 



This unique f(t) is said to be the representing function of the matrix mean 
M(A, B). So actually matrix means are in one to one correspondence with normal- 
ized operator monotone functions, the above characterization provides an order- 
isomorphism between them. Normalization means that /(l) = 1. For symmetric 
means, i.e. for means M(A,B) = M(B,A), we have f(t) = t/(l/t) which implies 
that /'(l) = 1/2. Operator monotone functions have strong continuity properties, 
namely all of them are analytic functions and can be analytically continued to the 
upper complex half-plane. This is the consequence of the integral characterization 
of an operator monotone function /(i), which is given over the interval (0,oo): 



where a is a real number, f3 > and fj, is a unique positive measure on (0, oo) such 
that 



Actually the interval (0, oo) may be changed to an arbitrary (a, 6), in this case the 
integral is transformed to this interval accordingly. These are the consequences of 
the theory of Loewner, an introduction to the theory can be found in Chapter 
V [5] . We will use this integral characterization at several points in the article in 
order to show that certain functions are analytic. 

In this article we are interested in finding all possible symmetric matrix means 
which are also geodesic midpoint operations on smooth manifolds. Or more gen- 
erally those matrix means that are arbitrary dividing points of geodesies. We will 
call such a matrix mean affine [32] : 

Definition 2 (Affine matrix mean) An affine matrix mean M : W 2 >-¥ W 
is a matrix mean which is also a point of an arc-length parametrized geodesic 
on a smooth manifold W 3 P(n, C) equipped with an affine connection V. I.e. 
M(A,B) = exp A (t\og A (B)) for a fixed t e (0,1) and for all A, B e P(n,C), B 
is assumed to be in the injectivity radius of the exponential map exp A (x) of the 
connection V given at the point A. The mapping \og A (x) is just the inverse of the 
exponential map at the point A g W. 








(G) 
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We can make some basic observations about affine matrix means. First of all 
note, that by (Jlj) we have that f(X) = M(I,X), so if M(A,B) is an affine matrix 
mean, then f(X) is some point of a geodesic connecting X and I. Also on a smooth 
manifold with an affine connection if we differentiate the exponential map exp p (X) 
at p, then we get d exp p = I p , where I p is the identity transformation of the tangent 
space at p [15] ■ Therefore if we differentiate its inverse, the logarithm map log p (g) 
we also get d log p = I p at p. So if we combine this with the chain rule we get that 
the differential of the mapping M(p, q) = exp p (t log p (q)) is dM(p, ■) = tl p . 

Now if we apply the above argument to an affine matrix mean M{A, B) we get 
the following result. 

Proposition 1 Let M(A,B) := exp A (t log^B)) be an affine matrix mean. Then 
/'(!) = *■ 

Proof Since P(n, C) is diffeomorphically embedded in H(n,C), therefore we can 
differentiate the map M(A, B) = exp^fi log^-B)) using the vector space structure 
of H(n, C), i.e. calculate the Frechet differential which we denote for an arbitrary 
differentiable function g by 

Dg[X][Y] = lim 9(X + Ys)-g{X) (?) 

at the matrix X in the direction of the matrix Y. So by ^ for all H € H{n, C) we 
have 

M(A,A + Hs) -M(A,A) 

s->0 S 

= lim A i/2 M(I, I + A~ 1 / 2 HA~ 1 / 2 s) - M(J, /) Al/2 = 

s->0 S 

= A 1 / 2 lim fil + A-^HA-^-fjl) 1/2 = A i/2 Df[I][A -V2 HA -l/2 ]A i/2. 

Since / is an operator monotone function on (0,oo), it admits an integral charac- 
terization (J3J) , so it can be analytically continued to the upper half- plane through 
the interval (0,oo). Therefore we may differentiate a power series representation 
of /, that uniformly converges on an open interval which contains 1, so then we 
get that Df[I][K] = Df[I][I]K = f'(l)K for all K G H(n,C)- Combining this with 
the above we get that 

^ M{A,A + Hs)-M(A,A) = ^ ^^-i/a^-i^ A H* = 

Since H was arbitrary this yields that t = /'(l), because dM(p, ■) = tl p and also 
the tangent space of P(n, C) at every point can be indentified by H(n,C). □ 

By the preceding proposition we shall focus on matrix means represented by 
operator monotone functions / on (0,oo) such that /'(l) £ (0, 1). We will use the 
notation ^P(t) to denote the set of all operator monotone functions / on (0, oo) 
such that f(x) > for all x G (0, oo) and /(l) = 1,/'(1) = t. We can find the 
minimal and maximal elements of ty(t) for all t G (0, 1) easily. 
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Lemma 1 For all f(x) £ we have 

((1 - 1) + tx' 1 ) 1 < f(x) < (1 -t) +tx. (8) 

Proof Since every operator monotone function is operator concave, see Chapter V 
[5], therefore we must have f(x) < (1 — t) + tx by concavity and the normalization 
conditions on elements of *P(t). Since the map x~ 1 is order reversing on positive 
matrices, we have that if f(x) £ ^P(t) then also /(x^ 1 )^ 1 € ^P(i). So again by 
concavity 

fix^y 1 < {i-t) + tx 

fix' 1 ) > ((l-^+tx)- 1 
f(x) > ((l-^+te" 1 ) 

□ 

Since ((1 — t) + tx~ x ) 1 and (1— t) +tx are operator monotone we see that the 
they are the minimal and maximal elements of ^J(t) respectively, and also they 
are the representing functions of the weighted harmonic and arithmetic means. 
This already gives us that the minimal and maximal affine matrix means are the 
weighted harmonic and arithmetic means respectively, so if M(A,B) is an affine 
matrix mean, then 

[(1 +tB~ 1 ] 1 < M(A,B) < (1 -t)A + tB. (9) 

In general by the previous Lemma [1] the above inequality is true for all M{A, B) 
matrix means with representing operator monotone function / for which we have 
/'(l) = t. In this sense *}3(i) characterizes weighted matrix means. If we take this 
as the definition of weighted matrix means, one can compare it with the definition 
of weighted matrix means given in |32j . 

Lemma 2 All f(x) g ^P(t) for t G (0, 1) has only one fixed point in (0, oo) which is 1 
and 1 is an attractive fixed point on (0,oo). 

Proof By the definition of ^{t) for all members f(x) of this set /(l) = 1, so 1 is 
indeed a fixed point. By the preceding Lemma [1] we have 

((1 -f) +faT 1 ) 1 < f{x) < (1 -t) +tx. 

Therefore for all x > 1 we have f(x) < x, i.e. f(x) has no fixed point in (l,oo). 
Similarly for all x £ (0, 1) we have 

x < ((l-^ + txr 1 ^ 1 <f(x), 

therefore f(x) cannot have a fixed point in (0, 1) as well. 

Now the attractivity of the fixed point follows from the fact that f{x) is mono- 
tonically increasing positive and concave on (0, oo) by operator monotonicity. Con- 
cavity implies that f"{x) < for all x £ (0, oo). Also its derivative /'(l) = t G (0, 1) 
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at the fixed point 1, so by Banach's fixed point theorem this fixed point is attrac- 
tive on (e, oo), where e < 1 is such that the derivative /'(e) = 1. On (0,e) the 
function (1 — t) + tx > f(x) > x so its subsequent iterates form an increasing se- 
quence of functions. I.e. if we start an iteration with xq £ (0,e), then after finitely 
many iterations by f(x), x n = f(x n —i) will be in the interval (e, 1). From there 
convergence to 1 follows again from Banach's fixed point theorem. □ 

In order to advance further in the understanding of affine matrix means, we 
should be able to grasp more geometrical structure related to the affinely connected 
manifolds corresponding to affine matrix means. In the next section we will study 
the general situation of affinely connected manifolds given with a geodesic dividing 
point operation. We will see that in this case we can reconstruct the exponential 
map and its inverse, the logarithm map from the geodesic dividing point operation. 



3 The reconstruction of the exponential map 

In this section we reconstruct the exponential map of an arbitrary affinely con- 
nected differentiable manifold based first on its midpoint map. Without loss of 
generality we fix a base point p as the starting point of the geodesies. The ba- 
sics of the exponential map of a manifold can be found for example in Chapter I. 
paragraph 6 [15] . 

Theorem 1 Let M be an affinely connected smooth manifold diffeomorphically embed- 
ded into a vector space V. Suppose that the midpoint map m(p,q) = exp p (l/2 log p (g)) 
is known in every normal neighborhood where the exponential map exp p (X) is a dif- 
feomorphism. Then in these normal neighborhoods the inverse of the exponential map 
logp(g) can be fully reconstructed from the midpoint map in the form 

, / \ r m(p,q) on -p 

logpW = lim — J . (10) 

where we use the notation m(p,q) on = m (p, m(p, g)°^ ,l_1 '^ . 

Proof We will use some basic properties of the differential of the exponential map 
to construct the inverse of it, the logarithm map. Since in small enough normal 
neighborhoods the exponential map is a diffeomorphism, it can be given as the 
inverse of the logarithm map log p (g). 

By the basic properties of the exponential map we have 



dexp p (Xt) 



<:)t 



exp(Xt)-p 
lim - = X, 



t=o 



where X e T P M. Here we used the fact that we have an embedding into a vector 
space. Suppose exp p (X) = q is in the normal neighborhood. We are going to 
provide the limit on the right hand side of the above equation. The limit clearly 
exists in the normal neighborhood so 

exp p (ja) -p = exp p (Jfjr) -P = m(p,q)° n ~p 



t—>Q t n— too — L. n—too — 
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Here we use the notation m(p,q) on = m yp,m(p,q) 0< - n 1 M. We are in a normal 
neighborhood so the exponential map has an inverse, the logarithm map, so the 
limit can be written as 

e XPp (Xt)~p m(p, g )°"- P . . 

X = } m A + = ^ 1 = lo Sp(9)- 

□ 

In the above assertion we used the midpoint map to reconstruct the exponen- 
tial map, but we can use arbitrary dividing point operation that yields a point, 
other then the ending points on the geodesic connecting two points in the normal 
neighborhood. This is summarized in the following proposition. 

Proposition 2 Let M be an affinely connected smooth manifold diffeomorphically em- 
bedded into a vector space V. In every normal neighborhood N let 7 a ^(t) denote the 
geodesic connecting a,b £ N with parametrization 7 Q ,{>(0) = a and j a b(l) = b. Suppose 
that the map m(a,b)t = 7 a ,fc(io) = ex P p (^o l°gp(?)) is known for a to £ (0,1) in every 
normal neighborhood N where the exponential map is a diffeomorphism and a,b 6 N . 
Then in these normal neighborhoods the logarithm map can be fully reconstructed as 



= hm m M%- P 



with the notation m(p,q)°™ = m (p,m(p,q) t ) . We also obtain the exponential 
map by inverting \ogp(q). 

We are going to use this construction in the next sections to characterize affine 
matrix means. 



4 The exponential map of affine matrix means 

Based on the idea of reconstruction given by Proposition[2]we are going to formally 
take the limits for matrix means in The following result will show that if 

a matrix mean is affine then the exponential map of the corresponding smooth 
manifold has a special structure. We will use similarly the notation M(A, B) on = 

as before in the previous section. 
Theorem 2 Let M(A, B) be a matrix mean with representing function f £ ^P(i) ■ Then 

J5S, M(A f B ( lyr A = ^' (.4-" 2 ^-" 2 ) ^ (11) 

where the limit exists and is uniform for all A,B£ P(n, C) and \ogj(x) is an operator 
monotone function which fulfills the functional equation 

lo g/ (/(x)) = /'(l)lo g/ ^) (12) 



on the interval (0,oo). 
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Proof We will prove the convergence to a continuous function logj(t) in a more 
general setting. The operator monotonicity in the matrix mean case will be a 
particularization. 

First of all note that by the repeated usage of ^ we can reduce the above 
problem to the right hand side of the following formula: 



07/ 



/'(I)" /'(I)" 

From now on we will explicitly use the notation g(x)° n = g (g(x)°^ n ~^J for arbi- 
trary function g(t) where this notation is straightforward. 

Due to the above formula it is enough to prove the assertion for a single operator 
monotone function f(x). By operator monotonicity of f(x) this is just the special 
case of the problem considered for arbitrary concave, analytic functions f(x) given 
in the following form 



for X G P(n, C). As every operator monotone function which maps (0, 00) to (0, 00), 
is analytic on (0, 00) and has an analytic continuation to the complex upper half- 
plane across (0, 00), we can consider the functional calculus for hermitian matrices 
in the above equations. Therefore we can further reduce the problem to the set 
of the positive reals by diagonalizing X and considering the convergence for every 
distinct diagonal element separately. For an extensive study on operator monotone 
functions one may refer to Chapter V in [5] • 

Without loss of generality we may shift the function f(x) by 1 so it is enough 
to show the assertion for 

n— *-oo g (0) n 

where g(x) = f(x + 1) — 1 and so g{x) on = f(x + l) on — 1. From now on we will 
be considering the shifted problem. At this point we must emphasize the fact that 
the function g must have as an attractive and only fixed point on the interval 
of interest (— l,oo). In the unshifted case this is equivalent to / having 1 as the 
only attractive fixed point on the interval (0,oo), which is the case by Lemma [2j 
So we can also assume that < g (0) < 1. The rest of the argument will be based 
on the claim that the above limit of analytic functions of the form g(x) on /g'(0) n 
is uniform Cauchy therefore the limit function exists and is continuous. 

First of all we have as the attractive and only fixed point of g, so for arbitrary 
x G (— l,oo) the sequence x n = g(x) on converges to 0. We have g(0) = and by 
the mean value theorem we have 

n 

x n = g(x)° n = g'{t n )g{x) o{n ~ 1) = Y\ 9'{k)x, 

i=l 

where t { G [0, ff(a;) o(i_1) ] if 1 > or i t e \^(x)°^ x \ oj if x < 0, since g is a 

concave function on (— l,oo). As x n — > for arbitrary x we have g'(ti) — > </(0). 
Now we have to obtain a suitable upper bound on 



gW" g(x) 



QUI. 



g'(0) n g'(0) r 



(14) 
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We argue as follows 



9{x)° r 



\g(xy n -g'(0) n - m g(xy 



S'(O)" g'(Q) m fl '(0)« 
„ \W=m + i9'(t t )-g'(0T- m \\IlZi9'(k)\ 
g'(0) n 

n f / \ rn t /_£ \ I 

n fi-i n#w 



i=m-|-l 



Now uniform convergence follows if \Y\iLi 9' (h) / 9' < 00 because then the tail 
rii^m+iS'^J/s'C) -> 1 so (|14[) can be arbitrarily small on any closed interval 
in ( — l,oo) by choosing a uniform m. By the continuity of g'(t) and x n — > we 
have g'(ti) — >■ </(0) and by assumption < g'{0) < 1, therefore there exists N and 
q such that for all i > iV we have < g'(tj) < g < 1. What follows here is that 
Bifi,^ < 00 such that |t jv I < Ki and \g"(U)\ < K 2 for all i > N. This yields the 
bound \ti\ < K\q l ~ N for all i > N . Considering the Taylor expansion of g around 
we get 

g'{k) _ </(0)+ 5 "(^ 



g'(0) 



ff'(0) 



for < t'i < ti. What follows from this is that 



n 



g'{U) 



fl'(0) 



n 



^1^2 i-jV 

1 + 7(oT g 



The infinite product on the right hand side converges because Yl'jLo g>(a) ^ con " 
verges hence [IlSi fl'(*i)/ff'(0)| < 00 f° r au ^ m the closed interval. 

At this point we can easily establish the convergence for normalized operator 
monotone functions because they are concave functions by Theorem V.2.5 in [5], 
so f"(t) < and they have only one fixed point which is 1. The fact that the limit 
is operator monotone function in this case follows from the operator monotonicity 
of the generating f(t). 

The functional equation (112[) is the consequence of the following: 



BIKJ ^ " n^oo f'(l) n n->oo /'(l) n 



□ 



Actually the above proof works for a larger class of functions then the family 
of normalized operator monotone functions. The limit in (|13p exists and it is a 
continuous function if the twicely differentiable function f(x) has 1 as the only 
attractive fixed point and the derivative — 1 < f'(x) < 1. The next example shows 
how to calculate the limit function explicitly. 
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Example 1 Consider the one-parameter family of functions 

f q (x) = [(l-t) + tx q ] 1/q 

for t G (0, 1). These are in ?p(t) if and only if q G [—1, 1], because for other values 
of q the function is not operator monotone, see exercise 4.5.11 in [6]. It is easy to 
see that 



n-l 



1/9 



[t n x q -t n + 1} 1/Q . 



t n x q +Y J tk v-- t ) 

k=0 

In this case we can easily calculate the limit function log/ f (x) by turning the 
limit into a derivative: 



, , .. (t n x q -t n + l) 1/q - 1 ,. (sx q - s + l) 1/( ? -1 
lo Si,f q ( x > = lim ^ — ^ = lim i 

n— >oo t s— >0 

l 



= ^-(sx q -s + l) 1 ' q 
as 

The limit functions indeed are operator monotone again if and only if q G [— 1, 1]. 
This family has a singularity at q = but it is easy to verify that it is a removable 
singularity, so in fact we have 

fo(x) = x l 
log/./oOO = i°g(z)> 

where log(a;) and x l are also well known to be operator monotone. Particularly x t 
as a representing function corresponds to the weighted geometric mean. 

Proposition 3 The limit function log/ (x) in Theorem fj| satisfies the following: 

(i) log/(a;) maps P(n,C) to H(n,C) injectively, 
(ii) 1 — x^ 1 < \ogj(x) < x — 1 for all x > 0, 

(Hi) If\ogj f(x) andlogj g (x) are the corresponding limit functions for f ', g G ^P(t) such 

that f(x) < g(x) for all x > 0, then log/ j(x) < log/ g (x) for all x > 0, 
(iv) log/(l) = and log'/(l) = 1. 

Proof (iii): Since f{x) < g{x) by monotonicity we have f{x) on < g(x) on . From this 
it follows that 

fjxT - 1 < qixT ^ 1 
" g'(l)» ' 

and the inequality is also preserved in the limit, 
(ii): By Lemma [1] we have 



({l-t) + tx ^ 1 < f(x) < (1 - t) + tx 



where on the left hand side we have the function and on the right hand 

side we have fi(x) from Example[TJ In Example[T]we calculated the corresponding 
limit functions, so these combined with the previous property (iii) proves property 
(ii). 

(i): By property (ii) it follows that log/ (a;) is nonconstant on (0,oo). Also 
log/ (x) is operator monotone there, so it is strictly concave, therefore injective 
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and real valued. This combined with the functional calculus for matrix functions 
proves the property. 

(iv): log J (l) = follows from (ii). Using this and (ii) again we have 

1-jl + h)- 1 < lo gj (l + fe)-lo gj (l) < (1 + /Q-1 
h h h 

Taking the limit h — ¥ we get derivatives on the left and right hand sides are 1, 
so also logj(l) = 1. □ 

Since \ogj(x) is operator monotone on (0, oo), it is also analytic there, so it has 
an analytic inverse exp/(a;) by Lagrange's inversion theorem, since its derivative is 
nonzero due to Proposition^ It is also easy to see that expj(0) = 1 and expj(O) = 
1. By these considerations we have just arrived at the following 



Proposition 4 Let f € <P(t). Then 



/(x)=ex P/ (/'(l)log 7 (x)), (15) 

where logj(x) is the unique solution of the functional equation (|12[) in the class of 
functions that are continuous, mvertible on (0, oo), vanish at 1 and have derivative 1 
at 1. 



Proof The first part of the assertion follows from the analyticity and invertibility of 
logj(a;) on (0, oo). For the second part note that if logj(ir) is an invertible solution 
of the functional equation (|12p and also log J (l) = and log'^l) = 1, then its 
inverse expj(:r) exists, exp 7 (0) = 1 and expj(0) = 1. Moreover 



lim 

n— »oo 

lim 

s-s-0 



\ogj(x). 



exp / (/ / (l)"lo g/ ( a; ))-exp / (0) 
/'(I)" 

exp J (slog 7 (x)) - exp 7 (0) _ d_ 
s ds 



exp J (slog J (o;)) 



\s=0 



□ 



The above propositions put some restrictions on the possible functions log/ (x) 
that can occur as limits in Theorem [2j Therefore we will use the notation £ to 
denote the set of operator monotone functions g(x) on (0,oo) such that g(l) = 
and <?'(1) = 1. By Proposition [3] it is clear, that for each / € ^P(t) we have 
a unique corresponding log^a;) in £. We will further say that for an arbitrary 
/ € ^P(t) the corresponding solution logj(s) of the functional equation (|12|) is the 
logarithm map corresponding to f(x), while its inverse exp/(a;) is the exponential 
map corresponding to f(x). 

In the following section we will go the other way around and see whether the 
function 

f(x) = ex P/ (t\ogj(x)) 



is in «P(i) for all log 7 6 £ and t € (0, 1). 
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5 Representing functions induced by logarithm maps 

In the previous section we established that for every / £ there exists a unique 
function log 7 £ £ such that it fulfills the functional equation (|12[) . In order to see 
whether an element \ogj £ £ also induces a representing function / £ ^P(t) with 
the generalized functional equation 



for all t £ (0, 1), we must extend our investigations into the upper complex half- 
plane H+ = {z £ C : Qz > 0}. 

First of all let us recall Nevanlinna's representation [5 of holomorphic functions 
/ : H + — > H + . By Nevanlinna's theorem each such / can uniquely be written as 



where a £ K, /3 > and v is a positive measure with support in (— oo, oo). It is well 
known that / can be extended to the lower half-plane HP = {z £ C : Qz < 0} as 
well by Schwarz reflection f(z) for all z £ HP. Therefore also if this extension is 
by analytic continuation over an interval (a, b) , then v vanishes on the interval [5] . 
Similarly if v vanishes on a real interval, then / is holomorphic on the interval as 
well and can be analyticly continued to the lower half-plane. 

Representation (jTHJ) will be useful for studying functions in £. For example 
Nevanlinna's representation yields that all / £ £ can be represented as 



where a £ R, /3 > and v is a positive measure with support in (— oo,0). This is 
due to the required holomorphicity of / on (0, oo) . Next let us find the maximal 
and minimal elements in £. 

Lemma 3 For all f(x) £ £ we have 



Proof Since every operator monotone function is operator concave, therefore we 
must have f(x) < x — 1 by concavity and the normalization conditions on elements 
of £. Since the map x^ 1 and — x is order reversing on hermitian matrices, we have 
that if f(x) £ £ then also —f(x~ 1 ) £ £. So again by concavity 



log 7 (/(x)) = t\og T (x) 




(16) 




(17) 



1 -x < f(x) <x-l. 



(18) 



fix- 1 ) <x-l 
f(x) >l-x-\ 



Clearly x — 1 and 1 — x 1 are also in £. □ 



At this point let us refer again to the functional equation (|12[) in the previous 
section. By the above considerations we can generalize (|12[) by analytic continua- 
tion. 



14 



Miklos Palfia 



Proposition 5 Let f G Then the function logj(a;) given in Theorem^ admits 

analytic continuation to H + and also to M~ across (0, oo) by relfection, moreover it 
fulfills the functional equation 

log I (/(*)) = / / (l)log I (z) (19) 

for all z G C\(-oo,0]. 

Proof Since analytic continuation of f(x) and logj(:r) can be performed using the 
integral characterizations ([5|) and (|17|) respectively, we end up with holomorphic 
functions living on C\(— oo,0]. Since \ogj(z) is a holomorphic function, it has a 
meromorphic inverse exp^(z). So we have 

F(z) = ex P/ (/'(l)lo gj (z)), 

a meromorphic function that is identical to f(z) everywhere on the domain (0, oo). 
Therefore by uniqueness of meromorphic and analytic continuation we must have 
F(z) = f(z) everywhere on the domain C\(— oo,0]. □ 

The above result tells us, that for a given logj G £ we should consider the 
generalized functional equation 

log J (/ t («))=tlog / (2) (20) 

to define a representing function ft{z) for all t G (0,1) corresponding to \ogj(z) 
which was itself obtained by analyitc continuation using representation (|17[) . The 
obvious question that arizes here is whether every logj(z) in £ has a corresponding 
ft G We need the following: 

Definition 3 (Radial convexity) Let SCCbe such that € S. We will say that 
S is radially convex if and only if for all z G S also tz 6 S for all t 6 [0, 1]. 

Proposition 6 Let log/ G £. Then log/ maps H + to a radially convex set in 

Proof Using Nevanlinna's characterization (|17[) we have that log/ is a convex com- 
bination of functions of the form 

, . \z + 1 
h,(z) = T —. 

If we have h\(z) = w for z £ H + , then after some calculations we get that 

z = X-(X +1) -2 
|A + w\ 

which means that 

A 2 + 1 



\X + w\ 2 ' 



In other words if h\{z) = w G H + , then for all s G (0, 1) there exists z s G H + such 
that h\(z s ) = sw. 

Now if we consider any convex combination of such functions h\(z), the re- 
sulting function will still have a radially convex image of H + . The reason for this 
is that if Xi G S, C H + where Si are radially convex sets, then sxj G Si for all 
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s £ (0, 1). Therefore if Si are the images of H + under the mappings Kih\.(z) for 
some K-i > and \, then the image S of H + under the function that we get as 
the sum of the functions Kih\ i (z), is radially convex, since every element of it can 
be written as a sum of some Xi £ Si. So we also have that the sum of sxi is in S 
too by the convexity of each Si. Therefore S must be radially convex. □ 

Theorem 3 Let log/ £ £. Then ft £ ^(t) for all t £ (0, 1) if and only if logj (z) has 
no branch point in 

Proof First of all since logj £ £, it follows that logj(x) is invertible on (0,oo) 
because it is nonconstant monotone increasing there, also it is invertible in a 
neighborhood of (0, oo) and its inverse expj(z) is holomorphic in that neighborhood 
and ft{x) £ (0, oo) and meromorphic in H + . 

Suppose that logj(z) has no branch point in H + . Then by the previous Propo- 
sition[6]it maps H + to a radially convex set. Since logj(z) has no branch point, it 
has a univalent holomorphic inverse expj(z), so 

ft(z) = ex Pl (t\ogj(z)) 

is a well defined holomorphic function on H + . Moreover ft{z) is real valued over 
(0,oo). Since the image log / (H+) of H+ under the map logj(z) is radially convex, 
we have that for any s £ log J (H + ) also ts £ log J (H + ). Therefore tlogj(EI + ) C 
logj(H+), so also exp J (tlog J (H+)) C M + . 

Now for the only if part suppose on the contrary that log/ has a branch point 
in H + . Then its inverse expj has a pole at the image of the branch point under logj 
which means that ft(z) is not holomorphic there, but this contradicts ft £ VP{t). 
□ 

What happens if log/ has a branch point in H + ? What can then be said about 
/t(*)? 

Proposition 7 Let logj £ £ be induced by an ft„ £ *P(£o) using Proposition^ Then 
ft £ ^P(t) for all0<t< to- 
Proof By Proposition 0] we have that there is no image of a branch point of log/ 
in the domain to \ogj(M + ) C H + , otherwise ft { z ) would have a singularity in H + . 
But since for all < t < to we have that tlogj(H + ) C to logj (H + ), therefore ft(z) 
is singularity free as well. □ 

Remark 1 In general one can assure that if for a given log/ £ £ with branch points 
t logj(H + ) avoids the image of the branch points (of log/) under logj in H + , then 

According to Theorem[3]we need to find members of £ without branch points. 
In other words we are looking for mappings that are univalent (schlicht) holomor- 
phic functions on HI" 1- mapping H + into itself. Such mappings are characterized by 
FitzGerald in the classical article [12] . 

Theorem 4 (FitzGerald) Suppose f(x) is a twice continuously differentiable, real- 
valued function with positive first derivative on (a,b). Suppose the origin is in (a,b) 
and /(0) = 0. A necessary and sufficient condition that f can be continued to be a 
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univalent analytic function of M + onto a subset of itself that is radially convex with 
respect to the origin is that the function 

V (x) - - - 



be conditionally positive definite, i.e. 

1 '•" s) Mz_M mdsdt > 

for all real continuous <f> having compact support in (a, b) and satisfying <j)(s)ds = 0, 
where Sifl^SW ts identified with t)'(s). 

Summarizing the results in the previous sections from the point of view of 
affine matrix means we arrive at the following 

Proposition 8 If a matrix mean M(A, B) is affine, then the exponential map and its 
inverse, the logarithm map are of the following forms 

I v\ 1/2 ( -1/2-v -l/2\ 1/2 

exp JX)=p< expj [p 1 Xp 1 )p 1 

, ( ( 21 ) 

log p (X)=p 1 / 2 log / (p- 1 / 2 Xp- 1 / 2 )p 1 / 2 

for p £ P(n,C), where expj(X) and \ogj(X) are analytic functions such that expj : 
H(n,C) M- P(n,C) andlogj(X) is its inverse andlog'j(I) = J,expj(0) = /,log/(7) = 
0,exp J (0) = I. 

Note that by Weierstrass's approximation theorem we also have 



1/2 

p 1 expj 



( -1/2 V -l/2\ 1/2 / -lv\ 

( p 1 Xp 1 I p 1 = p exp j [p X I 
V ' V ' (22) 

Vaiog, (p-^Xp-^)p^ =plo gl (p-'x) . 



In some cases, to ensure easier reading, similarly as in the above formulas, we 
will denote matrices with uppercase letters which are elements of some tangent 
space, while at the same time we will use lowercase letters for denoting matrices 
which are points of a differentiable manifold. 



6 Construction of an invariant affine connection 

Let us recall the classical symmetric space GL(n, C)/U(n, C), the cone of positive 
definite n x n matrices P(n, C) [9j . This is a Lie group and the K = U(n, C) 
isotropy group invariant inner product at the identity / is (U, V) = Tr {UV}. The 
tangent space, considering the Cartan decomposition of the Lie algebra, is the 
space of Hermitian matrices H(n, C) . The action of the isometry group GL(n, C) 
on this manifold is g[6) = gog* and acting with left translations we can transport 
the inner product to any point p on this manifold and we get the Riemannian 
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metric (U, V) = Tr \p 1 Up V}. The exponential map is just the ordinary matrix 
exponential at the identity. The left invariant affine connection is 

V Xp Y p = DY[p] [X p ] - i [Xpp-Wp + Ypp^Xp) , (23) 

here _DY[p]pf p ] denotes the Frechet-differential of Y at the point p in the direction 
Xp. A well known property of this metric is that the midpoint map of the space 
m (Pi?) = exp p (l /2log p (q)) is just the geometric mean of two positive matrices 

G(A,B) = A 1 ' 2 (A^^BA- 1 ' 2 ) 1/2 A 1 ' 2 . (24) 

The question that can be asked at this point is that are there other symmetric 
matrix means which correspond to symmetric spaces as midpoint maps on P(n, C)? 
Two other examples are known, these are the arithmetic mean (A + B)/2 and the 
harmonic mean 2(A~ 1 +B~ 1 )~ 1 . The symmetric spaces corresponding to these two 
means are Euclidean while the symmetric space corresponding to the geometric 
mean has nonpositive curvature. It has flat and negatively curved de Rham factors. 

At this point we begin with the characterization of means that correspond to 
affine symmetric spaces in general. What we know at this point is that the two 
functions, which are of each others inverse, logj(f) and expj(t) exist for all matrix 
means, as it was proved in Theorem [2j 

In [II] and [15] there is an extensive study of affine connections on manifolds. 
A well known fact is that the affine connection on a manifold can be reconstructed 
by differentiating the parallel transport: 

^ v .. /f( 7 )y 7(t) -y 7(0) 
v ** y > = j38 t ' 

where 7 (t) denotes an arbitrary smooth curve emanating from p in the direction 
X p = dj(t)/dt\t=o and rf(-y)Y denotes the parallel transport of the vector field 
Y along the curve 7 from j(t) to 7 (s), refer to [141115] . The above limit does not 
depend on the curve itself, only on its initial direction vector and it depends on 
the vector field Y in an open neighborhood of p. On affine symmetric spaces the 
parallel transport from one point to another along the connecting geodesic is given 
by the differential of the geodesic symmetries with a negative sign. The geodesic 
symmetry is given as 

S p (q) = exp p (-logp(?)). 

On affine symmetric spaces this map is an affine transformation so one can con- 
clude that 

zf ( 7 )y = - ag 7(i/ a )(«*«(")) , (25 ) 

" l t=0 

where j(t) is a geodesic connecting p = 7(0) and q = 7(1). 

We have already proved the following formulas for the exponential and loga- 
rithm maps at the end of the preceding section 



(26) 



expp(X) =p 1/2 exp / (p 1/2 Xp 1/2 jp 1/2 =pexp 7 (p 1 X 
log p (X) = P 1/2 log, (p^^Xp- 1 ' 2 ) p 1 ' 2 = plo gj (p^x) 
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The above identities already specify the geodesic symmetries with the notation 
Si(X) =exp J (-log / (X)) as 

S p (q) = ex Pp (- logp(«)) = P 1/2 Sj (p- 1/2 qp- 1/2 ) P 1/2 = P Sj (p-'q) . (27) 

Now we are in position to prove the following 

Theorem 5 Let P(n, C) be subset of an affine symmetric space with affine geodesic 
symmetries given as (|27[) . Then the invariant affine connection has the form 

V Xp Y p = DY[p] [X p ] - 1 (Xpp-Wp + Ypp^Xp) , (28) 

where k = S"(l)/2. 



Proof We are going to use (|25|) to obtain the connection (|28[) . We make the as- 
sumption that the geodesic symmetries are of the form (|27p . The functions exp p (X) 
and \og p (X) are of the form (|26|) . where expj(t) and logj(i) are analytic functions 
on a disk centered around and 1 respectively. We also have that log 7 (l) = 0, 
exp^(O) = 1 and furthermore 

dexpj(t) _ 
9t t=o ' 

First of all we have to differentiate the map S p (q) given in (|27p to obtain 
^1(7)^ = Tq-tpY, where j(t) is a geodesic connecting p = 7(0) and q = 7(1). 



dS p (expJYt)) 



dpSjip' 1 exp g (Yt)) 



dt 

pDSj [p- 1 ?] [ P - X Y 



t=o 9t 



t=o 



(29) 



We used the fact that <9exp g (Yi)/<9i|t = o = Y which is a consequence of exp'^O) = 1. 

Now we are going to differentiate the parallel transport as given by (|25|) to get 
back the connection. We use the holomorphic functional calculus to express the 
Frechet-differential in (1291) as 

DS I \X}\U} = ^- I SAz)[zI - X^UlzI - X^dz. 

2m J 9 

It also easy to see that 7?S/[7][7] = <S/(1) = — 1, so we may express the limit (|25|) 
by the following differential 

a^(t/2)DS! ^mr 1 ^)] [y{t/2)-% {t) ] 



Vy (0)^7(0) ; qI 



t=o 
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we massage this further by using the holomorphic functional calculus 



S T (z)[zI - 7 (t/2)- 1 7W]" 1 7(i/2)" 1 V 7(t) x 



5 7- 7 (t/2)- 1 7(i)]" 1 & 



1 , 



-sY(0)7(0)~%(o)£>S'j[J][J]- 



+ [zl - ir^iOy'Y^izI - J]- 1 I 7 (0)- 1 7'(0)[ 2 / - !}-' + 



[zl-l]- 1 } 



7(0)- 1 i 7 '(0)7(0)- 1 y 7(0) + 7(0)- 1 W[ 7 (0)][ 7 '(0)] 



dz 



by using the fact that DSj [I] [I] and [zl — I] 1 commutes with every matrix we 
get 

. -™y ( o) 7 (o,-V 7(or 

- 7 <°>^t^>rV(oworv, (or 

~ 7(0) 2S [^ 7<orl7 ' (oh(orIr,<0 ^ 

at this point we use the integral representation 

to further simplify the above. 

V 7'(0)^r(0) = -^J^ [7'(0) 7 (0)-V 7(0) + 5^(0)7(0)- V(0)] - 

- ^7'(0) 7 (0)-V 7(0) - ^1 [-7'(0) T (0)-V 7(o) + 2W[ T (0)][ 7 '(0)]| = 



^(1)W[7(0)][ 7 '(0)] - ^ [ 7 '(0)7(0)-V 7(0) +y 7(0) 7(0)-V(0)l 



So we have that k = S"(l)/2. □ 

The above clearly tells us that all symmetric spaces occuring in such a way that 
their geodesic division maps are matrix means, have invariant affine connections 
in the form (|28p . We are going to study these connections as re being a parameter. 
We will find out later for which values of re are these spaces symmetric. Also for 
arbitrary real re (|28|) defines an affine connection with corresponding exponential 
and logarithm map which are of the form (|26|) as we will see later. We will also 
determine if these connections are metric or not. 
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7 Properties of these affine connections 

In this section we study the connections 

V X Y P = DY[p] [X p ] - % (Xpp-% + (30) 

for p G P(n, C) and vector fields X p , Y p e H(n, C) on the smooth manifold P(n, C) 
with tangent bundle H(n,C). It is easy to see that indeed these connections are 
affine and analytic for real k. 

We can fix a coordinate frame by taking the basis E^ G H(n, C) , where i in- 
dices over the set of distinct hermitian matrices which have zero entries, excluding 
exactly the entry = 1 and its transpose [-Ej];fc = 1. If we equip H(n,C) with 

the inner product (X, Y) = Tr{XY}, then the E^ form an orthonormal basis of 
H(n, C). The dimension of H(n, C) is n(n+l)/2 such as the dimension of the smooth 
manifold P(n, C) . In this coordinate frame the Christoffel symbols are given as 

r?jE k = ~ (e.p^Ej + Ejp^E^j , (31) 

where we used the Einstein summation convention for repeated covariant and 
contravariant indices. Given an arbitrary connection V the geodesic equations 
corresponding to it are given as 

V-Ht)7M = (32) 

with given initial conditions 7(0) and 7(0), for all t £ [0,T). I.e. the curve 7 must 
be parallel along itself. 



Theorem 6 The geodesic equations corresponding to the affine connections (|30[) are 

7 = K77 (33) 

The solutions of these equations with initial conditions 7(0) = p, 7(0) = X are the 
folowing one parameter family of functions 

V2 „™ (r,- 1 / 2 Y n - 1 / 2 A r, 1 / 2 



7 (t) = exp p (Xi) = p L ^ ex P/ [p-^Xp-^t) p 1 ^ (34) 
where 

exp^^-f + (35) 
lexp(Aj else. 

Proof For the connections (|30|) it is easy to see that the corresponding V^,( t )7(t) = 
geodesic equations are ([33]) . 

Let us first consider the case when 7(0) = p = I = 7(0) = X. Then it is enough 
to solve the equation (|33[) for real numbers. Therefore the equation takes the form 

expj'(t) = k exp'j(t) 2 expj(i) -1 . (36) 

If we transform the equation to the inverse function of expj(i) which will be the 
logarithm map logj(f), then we get a separable first order differential equation of 
the form 

iog , ;w = -^io g ;(^- 1 . 



Classification of affine matrix means 



21 



Solving the above we get the logarithm map as 



lo gj (X) 



1-K 

log(X) else 



if rc^l, 



From this by inverting the above function we get the assertion for real numbers. 
Now we check by substitution into (|33|) that the curve 

^/ 2 o^ L- 1 /2y„-i/2,L 1/2 



•y(t) = p expj I p ' Xp tjp 



is also a solution of the equations (|33[) . since the function expj is analytic. Indeed 

■ U\ V -1/2 ' I —1/2 V -1/2 A 1/2 1/2 / ( -1/2 ir -1/2 A -1/2 v 

7(t) = Ap ' expj [p ' Xp ' t\p ' =p' expj Ip ' Xp t\p ' X 



X 



j(t) = Xp- 1 / 2 exp? (p-^Xp-^h) p 
and after substitution we get 

v -1/2 11 I -1/2 v -1/2A -1/2 V -1/2 / / -1/2 v- -1/2, \ 

Xp ' exp ? ( p ' Xp t\p ' X = nXp ' exp ? ( p ' Xp 'tlx 



exp, 



( -1/2 v -1/2 A ' ( -1/2 v -1/2 A -1 

( p ' Xp ' t 1 exp j ( p 1 Xp ' t\p 



/2 



which is fulfilled since 



// / -1/2 v -1/2. 
expj (p ' Xp 1 t 



n expj 



-l/2 v -1/2,1 
' Ap ' t 1 expj 



- 1 / 2 xp- 1 / 2 f 



holds by the functional calculus for exp 7 and its derivatives and (|36[) . □ 

Corollary 1 The exponential and logarithm map for the affine connections (|30|) are 
given in the form 



e Wp (X)=p^e XPl (p-^Xp-^)p^ 
log p (X)=p 1 / 2 log / (p- 1 / 2 Xp- 1 /2) p i/2 i 



(37) 



where 



expj(Jf) 



log 7 (X) 



exp(X) eZse, 
log(X) eZse. 



(38) 



The affine matrix means which induce these affinely connected manifolds are 
M t (X,Y) =ex Pjf (tlog x (r)) = 



X l/2 
X l/2 



(1 -t)I + t(X- 1/2 YX 



X~ 1/2 YX~ 



-1/2 



X 1 / 2 ifn^l, (39) 



1/2) X l/2 



else, 



if K £ [0,2], /or ot/ier values of k the functions (|39[) are not matrix means. 
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Proof The first part of the assertion is clear, the second part follows from the fact 
that (|39[) are matrix means if and only if re G [0, 2], refer to Example [1] □ 

Remark 2 The one-parameter family of affine matrix means (|39[) seems to have a 
singularity at re = 1, however it is known that the singularity is removable and 
indeed as re — > 1 we get the matrix geometric mean as the limit. This phenomenon 
has already been investigated in [26 . In that paper the same one-parameter family 
of matrix means were considered under the name of matrix power means. 

If re = we get back the arithmetic mean as the midpoint operation, and the 
weighted arithmetic mean 

At(A, B) = (1 — t)A + tB (40) 

is the geodesic line connecting A and B with respect to the metric (X, Y) = 
Tr {X*Y}. If re = 2 we get back the harmonic mean as the midpoint operation, 
and the weighted harmonic mean 

H t (A,B)=((l-t)A- 1 +tB- 1 '\ 1 (41) 

is also a geodesic with respect to the metric (X, Y) p = Tr {p~ 2 Xp~ 2 Y} . We have 
already mentioned that the second metric is isometric to the first one, so it is also 
Euclidean. 

In the case when re = 1 the midpoint is the geometric mean and the geodesies 
are given by the weighted geometric mean 

G t (A, B) = A 1/2 (a- 1/2 BA- 1/2 ^ ' A 1 ' 2 . (42) 

The corresponding Riemannian metric is {X, Y) = Tr {p~ 1 Xp~ 1 Y}. This mani- 
fold, which is the symmetric space GL(n, C) / U(n, C) , is nonpositively curved while 
the other two has zero curvature. 

In the paper |26j matrix power means P s (wi, . . . , w^; A%, . . . , A^) are defined as 
the unique positive definite solution of the equations 

A' 

X = Y J W i G s (X,A i ) (43) 

i=l 

where s £ [—1, l],iOj > 0, Y^ii=i w i = 1 ana * e P(n,C). Existence and uniqueness 
of the solutions follow from the fact that the function 

k 

f(X) = J2^G s (X,A i ) 
i=l 

is a strict contraction for s£[-l,l],s^0 with repsect to Thompson's part metric 
[26] . In the case k = 2 we get back the affine matrix means (|39[) with s = re — 1 
and t = u>2- 

Corollary 2 With the identification s = re — 1, the two-variable matrix power means 
Ps{w i, 102 ; A\, A2) are geodesic lines, with 102 being the arc-length parameter, of the 
affinely connected spaces with affine connections (|30|) • 
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The arithmetic (|40[) , harmonic (|41[) and geometric (|42p means have nice chrac- 
terizations and extensions to several variables as the center of mass or Karcher 
mean of the corresponding manifolds [31,32,7,26 . I.e. 

n 

A(wi, . . . ,w k - Ax, . . . ,A k ) = argmin V] Wid 2 (X, A;), (44) 

where d(-, ■) is a Riemannian metric given as 

d 2 (X,Y) = (\og x (Y),\og x (Y)) x (45) 

where (■. ■) x is one of the corresponding metrics given above for the arithmetic 
(|40p . harmonic (|4ip and geometric (|42|) means, and \og x (Y) are the correspond- 
ing logarithm maps (|37|) (for « = 0,2,1 respectively). It is well known that in 
geodesically convex neighborhoods on a Riemannian manifold (|44[) has a unique 
solution [21,31]. The solution can be expressed by taking the gradient of the cost 
function on the right hand side of (|44|) [2"T] and then one arrives at 

n 

^2wilog x (Ai)=0. (46) 

i=l 

The unique solution of this equation can be expressed in closed form in the case 
of the arithmetic and harmonic means, since the corresponding manifolds are Eu- 
clidean. The solutions are just the multivariable weighted arithmetic >,-_-, lUj-Aj 

and harmonic means fj^i=i w i j! ^7 1 ^j E2- These functions are monotone in their 
variables with respect to the positive definite order and have some other desirable 
properties [26]. These two cases are well known and of less interest, however the 
same situation is of much more interest in the case of the geometric mean. In this 
case the unique solution of the minimization problem (|44[) cannot be expressed 
easily in closed form since the corresponding Riemannian manifold is no longer 
flat. The corresponding equation for the gradient (|46p is given in the form 

n 

J2^g(X 1 / 2 Ar 1 X 1 / 2 ) = (47) 

z=l 

and usually this equation is called the Karcher equation [26 and the corresponding 
unique positive definite solution (|44[) the Karcher mean. Several properties of this 
mean were open problems, for example its monotonicity with repsect to the positive 
definite order, however this and other key properties of the mean were proved by 
using different techniques [25 , 26||8] . One of the techniques given in [26] is based on 
the matrix power means P s (wi , . . . , Wk ;Ai,..., A^). These means are given as the 
unique positive definite solutions of (|43p . The following result provides a geometric 
characterization of matrix power means. 

Proposition 9 The matrix power means P s (w\, . . . , w/.; A±, . . . , A/.) for s 6 [—1, 1] 
are the unique positive definite solutions of the Karcher equations 

n 

J2mlog x (A t ) = (48) 

i=l 

where \og x (Aj) are the logarithm maps (|3T|) corresponding to the affine family (|30[) 
with parameter identification s = k — 1. 
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Proof The defining equation (|43l) of matrix power means with s ^ is equivalent 
to 

k 
i=l 

n 

^^x 1 / 2 [(x- 1 ' 2 ^- 112 ) 8 - i] x 1 ' 2 

i=l 

(x^^AX- 1 ^ 2 ) 3 -i 

J^WtX 1 / 2 ^ 1 X 1 / 2 

•i=l 

which is by (|37[) equivalent to 

n 

8=1 

The case s = is just the case (|47p . □ 

Remark 3 By the continuity of fixed points of pointwise continuous families of strict 
contractions [26] , the unique solution of (|48[) varies continuously with resppect to 
the parameter s. The singularity at s = is known to be removable and the limit 
is just the Karcher mean [26] . 

Now since the Karcher equations (|48|) admit unique positive definite solutions 
the obvious question arises whether there are Riemannian metrics corresponding to 
other values of k? Also for what other values of k is the manifold P(n, C) with affine 
connection (|30p a symmetric space? The full solution of these questions requires 
the study of the curvature tensors and holonomy groups which is postponed to the 
last section. At this point we prove some other results which gets us closer to this 
metrizability problem of the affine connections (|30[) . First of all we compute the 
parallel transport over a geodesic connecting an arbitrary point and the identity. 
The parallel transport of a vector Y 7 (o) given in the tangent space at 7(0) with 
respect to the connection V along the curve 7(t) is defined to be the vector field 
V 7 ( t ) which is the solution of the ODE 

V 7(t) F = 0. 

Proposition 10 Let c(t) be a geodesic with repsect to the connection (|30p and c(0) = 
I,c{l) = p. Then the unique solution ofX/^^Y = with respect to the connection 
(|30p and the initial condition Y c ^ = Yq is the vector field 

y(t) = c(t)*Y c(t)*. (49) 
Proof We have to integrate the equation V ' c / (ty^c(t) = 0. This is equivalent to 

DY[c(t)][c'(t)] - 1 ( C '(t) c (t)-V c(t) + r c(t)C (t)-V(t)) = 0. 

We are looking for the solution Y c ^ = Y(t) in the form 

Y(t) = f(c(t))Y f(c(t)), 



= 
= 

= 0, 
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for some analytic function f(x). We have for the Frechet-differential 

Now substituting into the equation of the parallel transport above, we get 

Y /(c(i)) + /(c(i))Y ^P = I (c'Wc(rV cW + K c(t) c( t )-y (t)) . 

Since c(t) = exp^(i logj(p)), it has a power series expansion, as has f(x), so we 
have by commutativity that 

^{t)c{t)- X f{c{t)) = = D/[c(t)][c'(t)] = /'(c)c'(i). 

Since everything on the left and right hand side commutes with one another, we 
arrive at the following separable differential equation 

which has its solution in the form /(c) = c K / 2 . □ 

On a Riemannian manifold the length of vectors are left invariant by parallel 
transports with respect to the Levi-Civita connection due to the Fundamental 
Theorem of Riemannian geometry [22]. It is easy to check that the connections 
(|30p are symmetric and torsion free so any of them can possibly be a Levi-Civita 
connection of a Riemannian manifold. So by the above proposition we should look 
for the Riemannian metrics in the form 



-«/2y„-*/2 -k/2 v -k/2 



(50) 



for some positive definite bilinear forms (•, ■) given on the tangent space at /. In 
the next section we prove that all affine matrix means are actually matrix power 
means, i.e. we do not have to look for other connections than (|30[) . 



8 The classification of affine matrix means 

Due to Proposition [S] we have the exponential and logarithm map of affine matrix 
means in the form 

fv\ 1/2 I —l/2v -l/2\ 1/2 

exp (X)=p / expj (p 1 Xp ' )p> 

) \ ( 51 ) 

\o gp (X)=p^lo gl (p~^Xp-^)p^ 

for p £ P(n,C), where exp 7 (X) and logj(X) are analytic functions. The func- 
tion expj : H(n, C) h-> P(n, C) and log/(X) is its inverse, log/(J) = /,expj(0) = 
I,logj(I) = 0,expj(0) = /. Suppose that (|5ip represent the exponential and loga- 
rithm map of an affinely connected manifold. Then the analytic function expj(t) 
is the solution of some geodesic equations 

exp'jr'(t) + r (expj(t),exp^(t),exp J (f)) = 0, 
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where _T(-,-,-) : H(n,C) x H(n,C) x Pin, C) M- H(n,C) is a smooth function in 
all variables and linear in the first two, representing the Christoffel symbols of an 
affine connection. By Propostion 15 and Corollary 16 of Chapter 6 in [33J we know 
that connections which have the same torsion and geodesies are identical and for an 
arbitrary connection there is a unique connection with vanishing torsion and with 
the same geodesies. If we have an affine connection with non-symmetric Christoffel 

symbols r^ k , it has the same geodesies as its symmetric part Jk 2 k3 , so without 
loss of generality we can assume in our case that all connections are symmetric, 
so we will be considering mappings _T(-, •, ■) which are symmetric in their first two 
arguments. 

Proposition 11 Suppose that _T(-, •, •), expj(-), exp p (-) are functions given with the 
above properties. Then 

r{x,x, P )= P 1 ' 2 r(p- 1 ' 2 x p - 1 ' 2 ,p- 1 / 2 Xp- 1 ' 2 ,i) P 1 ' 2 (52) 

for p e P(n, C) and X € H{n, C) . 

Proof Let 7(t) = expj fp — 1//2 Xp^ 1 ^ 2 ^ . Since expj is an analytic function we have 

■ / ,\ -1/2 v -1/2 / / -1/2 V -I/2A 

j(t) = p ' Xp ' expj [p ' Xp ' t\ 

■■/,\ -1/2^ -1/2 // / -l/2v -I/2A -1/2 V -1/2 

j(t) =p ' Xp ' exp/ \p ' Xp t\p ' Xp ' 

and other formulas hold for ^(t) and 7(i) similarly to the second part of the proof 
of Theorem [6] By the geodesic equations we have 

7(*) = --T(7(*),7(*),7(*)) 
Xp-^exp'/ (p- 1/2 Xp- 1/2 iV 1/2 X = - P 1/2 r (p-^Xp-^exp'j (p-V 2 x 

v -l/2,\ -l/2 v -1/2 / / -1/2 V -l/2,\ / -l/2v -l/2,\\ 1/2 

Xp 1 t\ ,p 1 Xp ' expj ( p ' Xp ' 1 1 , exp 7 ( p ' Xp ' 1 1 I p ' . 
If we consider the geodesic equations for 7(t) = exp J9 (Xt) we get 

v -1/2 III -1/2 v -l/2,\ -1/2 v ( v -1/2 ' ( -1/2 v -l/2,\ 1/2 

Xp ' expj (p ' Xp t\p ' X = ~T [Xp ' expj \p 1 Xp ' t\p ' , 

1/2 ' / -1/2 V -l/2,\ -l/2v 1/2 / -1/2 v -l/2,\ l/2\ 

p ' expj ( p ' Xp t\p ' X,p ' expj ( p ' Xp ' t\ p ' I . 

The left hand sides of the two equations above are the same so as the right hand 
sides. Taking t = and that expj(O) = 7, exp 7 (0) = I we get for allp € P(n, C),X 6 
H(n, C) that 

p 1 ' 2 r{p- 1 ' 2 xp- 1 / 2 ,p- 1 / 2 x p - 1 ' 2 j) p 1 ' 2 = 

= T(X,X,p), 

which proves the assertion. □ 

By the above result we have just reduced the problem of characterizing r (X, X, p) 
to the characterzation of r {X, X, I). Now we will show that r (X, X,p) is invariant 
under similarity transformations. 
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Proposition 12 For all p £ P(n, C) and X € H{n, C) and invertible S we have 

r (sxs- 1 ,sxs-\s P s- 1 ^ = sr(x,x, P )s~ 1 . (53) 

Proof We have by the geodesic equations 

X 2 exp"(Xt) = -r(Xexp' I (Xt),Xexp' I (Xt),exp I (Xt)) 
SX 2 exp"(Xt)S~ 1 = -Sr (X ex.p'j(Xt),X ex.p'jiX^^ex.pjiXt)) S' 1 . 

Similarly if we consider the geodesic equations for the curve j(t) = exp 7 (SXS^ 1 ^ 
we get 

SX 2 S~ 1 expiiSXS'H) = -r (SXS^ 1 exp' 7 (SXS _1 t), SXS^ 1 expiiSXS'h), 
ex.p I (SXS~ 1 t)\ 
SX 2 ex.p"(Xt)S~ 1 = -r ^SXex.p' I (Xt)S~ 1 ,SXexp' I (Xt)S~ 1 , 
Sexpj(Xt)S^ . 

Again since the above two equations are identical we get the assertion. □ 

By the above proposition we have for hermitian X that 

r(X,X,I) = Ur(D,D,I)U*, (54) 

for some diagonal D and unitary U, so it is enough to characterize F (X, X, I) for 
diagonal X. 

Theorem 7 Let D be diagonal with real coefficients. Then 

r(D,D,I) = -cD 2 , (55) 

for some real valued constant c. 

Proof First we will show that r (J, /, I) = cl for some real constant c. Consider 
the case when -y(t) = exp 7 (A/t) for some real A. Then by the geodesic equations 
for j(t) we have 

A 2 exp"(A7V) = -r(Xexp I (XIt),Xex.p' I (XIt),exp I (XIt)) . 

By linearity of P(-, •, •) in the first two variables, this is equivalent to 

A 2 exp"(A/i) = -A 2 r (expj(AJt),expj(A/t),expj(AJt)) • 
Letting ( = 0we get 

a = -r (1,1,1), 

where c = exp"(0) is a real number, since exp 7 : H(n,C) n> P(n, C) is an analytic 
function with real coefficients in its Taylor series. 

The next step is to show that for a projection P = P 2 = P* we have r (P, P, I) = 
—cP. Consider again ^(t) = exp/(Pt). Then the geodesic equations read 

P 2 exp"(Pt) = -T (Pexp' I (Pt),Pex.p' I (Pt),exp I (Pt)) . 
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Since P 2 = P and again letting t = we get 

cP = -r(P,P,i), 

where c is trivially the same constant as determined above for r (1,1,1). Now 
suppose that we have two mutually orthogonal projections Pi , P2 such that P1P2 = 
0. Then we have for the projection Pi + P2 using linearity of P(-, •, •) in the first 
two variables that 

P (Pi, Pi, /) + P (P 2 , P 2 , 1) = -c(P a + P 2 ) = P (Pi + P 2 , Pi + P 2 , J) = 
= P (Pi , Pi , /) + P (Pi , P 2 , 1) + r (P 2 , Pi , 7) + P (P 2 , P 2 , 1) , 

which yields that for mutually orthogonal projections Pi,P 2 we get the orthogo- 
nality relation 

P(Pi,P 2 ,J) =0. 

Finally since a diagonal D can be written as D = ^\ AjPj for mutually orthog- 
onal projections Pi, we have 



p (d, D,i) = r XiPi, *i p i> 1) 
= {p i ,p i ,i) = -J2 A-cPi = 



-cD 2 , 



which is what needed to be shown. □ 



The above three theorems with the other preceeding results presented here, 
lead us to the concluding 



Theorem 8 All affine matrix means Mt(X,Y) are of the form 
X 1 ' 2 



M(X,Y) 



X 1/2 ifn^l, 

(56) 



(l-^I + t^X-^YX- 1 / 2 

X 3.,* , x -l/2y X -l/2\* x l/2 ifK=l, 

where < k < 2. The symmetric affine connections corresponding to these means are 
V Xp Y p = DY[p] [X p ] - ^ (Xpp-Wp + Ypp^Xp) . (57) 

Proof By Proposition II II 1121 and Theorem[7]we have that the functions P(-, •, •) : 
H(n, C) x H(n, C) x P(n, C) H(n, C) representing the Christoffel symbols are of 
the form 

P{X,X,p) = -cXp^X. (58) 

This formula determines the functions that are the symmetric parts of the possible 
connections, and these connections have geodesies determined by Theorem[6]in the 
form (|56p . Again by Propostion 15 and Corollary 16 of Chapter 6 in [33] we know 
that connections which have the same torsion and geodesies are identical and 
for an arbitrary connection there is a unique connection with vanishing torsion 
and with the same geodesies. So in other words since the connections (|5T|) are 
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symmetric, affine and have the same geodesies, therefore they give the sought 
symmetric connections for each re if we choose c = re. 

The corresponding geodesies are given in (|39p . and these are matrix means if 
and only if re g [0,2], since the representing functions f(t) in turn out to be 
operator monotone only in these cases due to Example [TJ □ 

The above result gives us the complete classification of affine matrix means. So 
now we can concetrate only on the connections (|57p . In the next section we solve 
the metrization problem of these connections. 



9 The holonomy groups and metrizability of the affine family 

Let W be a smooth connected manifold with an affine connection V. The holonomy 
group "Hp(V) of the connection V at point p £ W is defined to be the set of all 
linear automorphisms of the tangent space T P W at p induced by parallel transports 
along p based closed rectificable curves. If W is simply connected then H P (V) is 
known to be a Lie subgroup of End(T p W / ) [27 . In case of non-simply connectedness 
the restricted holonomy group 'Hp(V) is defined as the normal subgroup of H P (V) 
which is induced by closed rectificable curves homotopic to zero, see Chapter II 
Section 4 in [52] for more detailed information. Let f)p(V) and f)p(V) denote the Lie 
algebra of "H P (V) and H P (X7) respectively. The holonomy group H P (V) is known 
to be an invariant of the connected manifold W, since H P (V) is conjugate to every 
other H q (V) by parallel transports. 

Now suppose that the connection V is real analytic. Then by Theorem 10.8 
of Chapter II and Theorem 9.2 of Chapter III in [22], f)p(V) is generated by the 
successive covariant differentials V r R, r = 0, 1,2, . . . at the point p where R(X, Y) 
denotes the curvature endomorphism of the connection V. This is a version of 
Ambrose-Singer's theorem of Kobayashi-Nomizu. The curvature tensor R is defined 
as 

R(X, Y)Z = Vx^yZ - V Y V X Z - V [x ,y]Z 
or expressed in local coordinate system with the Christoffel symbols 7^ fe as 

f)pi a pi 

R i _ Ij kj r i r m r i r m /r Q \ 

K 3 kl ~ ~Q^k faT ' km '3 ~ lm k 3 ' ' ■ > 

Suppose now that the connection V is torsion-free, i.e. 

V x y-V y X-[X,F]=0 (60) 

for all vector fields X, Y or equivalently _T* fe = rjy everywhere. Then V is the 
Levi-Civita connection of a Riemannian metric if and only if the corresponding 
holonomy group H P (S7) is a compact Lie group. More generally there exists a 
non-degenerate V invariant bilinear form (•, •) if and only if Hp(V) leaves (■, ■) 
invariant. 

In [27] all possible irreducible holonomy groups of torsion-free affine connec- 
tions are classified, so in principle we know what kind of groups can occur, at least 
in the reducible case. Again we are interested in the connections 



V Xp Y P = DY[p][X p ] - | (Xpp-% + Ypp^Xp) 



(61) 
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These connections are real analytic, torsion-free and the corresponding manifold 
P(n, C) is analytic simply connected. So to answer the question of metrizability 
we have to determine the holonomy groups H P (V). 
In our case it turns out that 



R U E i = (f - t) p [ p_1 ^' [p~ 1e *>p~ 1e i 



(62) 



where the form the standard basis of the vector space of H(n, C) and [■, ■] is the 
commutator. Note that the tangent space is H(n,C), so the left hand sides are in 
H(n, C). In order to determine which of these manifolds are symmetric spaces it is 
sufficient to calculate the covariant differential Rjf-i. m , since it vanishes everywhere 
if and only if the underlying manifold is a symmetric space [15] . Given the basis 
Ei for H(n, C) we have the identities 

R(X,Y;A 1 ,...,A r )Z = £ R] km . h _ lr E l X 3 Y k Z m A l 1 1 ---A l { 

i,j,k,l-i ,...,/ r ,m 
fi 

p! _ u p» I pi r>s 

f^jkmih ,...,l r +i Q x i r+1 rL jkm;l 1 ,....l T T 1 sl r + 1 rt jkm\l 1 .. . .,l r 

~ r'l IJt l r+1 Rjkm;h,...,s,...,lr 
A" 

where indices after ; denote covariant differentiation. Now we prove an analogue 
of Lemma 1 given in the proof of Theorem 9.2 of Chapter III [22] , 

Theorem 9 Let the smooth connected manifold P(n, C) be equipped with real analytic 
connection V and curvature tensor given by (|62[) with s£R. Then 



R(X, Y; A lt . . .,A r )Z = (1 - k) t D(R(X, Y)Z)\p)[A x , ...,A r ) (63) 

where D(R(X, Y)Z) [p] [A\, . . . , A r ] denotes the r-th Frechet differential of the map 
R(X,Y)Z at the point p £ P(n, C) in the directions Ai G H(n,C). 

Proof The proof is based on writing R(X, Y)Z and its subsequent covariant differ- 
entials in essentially two equivalent ways. First of all note that 

Ap-i = D{x- X )\p][E t ] = -p-'E.p- 1 , (64) 

so the differential operator is equivalent to Frechet differentiation at p in the 
direction of E4, also 



R{X,Y)Z=^-^y[p-'z\p- 1 X,p-'Y] 

^pj [Z [p- x X,p- x Y\ + \Yp-\Xp~ 1 ] Z) (65) 
[zp-\ \xp-\Yp- 
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Using index-less notation and the linearity of R(X, Y; A\, . . . , A r )Z we have 
R(X, Y;A X ,..., A r+1 )Z = V Ar+1 (R(X, Y; A\, . .., A r )Z) 

- ~ Ur+iP^RiX, Y; Ai, . . . , A r )Z + R(X, Y ; A u . . . , A r )Zp~ X A r+r 

- RiAr+^X + Xp~ 1 A r +i,Y; A u . . .,A r )Z 

- R(X, Ar+ip-W + Yp- 1 A r+1 -A 1 , ...,A r )Z (66) 

- R(X, Y\Ai,..., Ar^Ar+ip^Z + Zp^Ar+i) 

- J2 R ( x > Y;Ai,..., Ar+jp^Ai + Afp^Ar+i, A r )Z 1 . 
i=l J 

Again the first term in the above equation is equivalent to 

(A r+1 ) S ^R(X,Y; A x , . . .,A r )Z = D(R(X,Y; At,. . . , A r )Z) [p] [A r+1 ]. (67) 

Claim R(X, Y; A\, . . . , A r )Z is the linear combination of terms pS, where S is some 
word which is a product of the terms p~~ 1 X,p~ 1 Y,p~ 1 Ai, . . . ,p~ 1 A r of the first 
order. 

Proof (of the claim) We prove by induction. For r = it clearly holds by the first 
equality in (|65|) . Suppose that it holds for some r. Then by (I66[) it is easy to see 
that it holds for r + 1, due to (|64[) . the linearity of R(X, Y; A±, . . . , A r )Z and the 
product rule of Frechet differentiation. 

By the claim R(X, Y; A±, . . . , A r )Z is the linear combination of terms pS, therefore 
by linearity, (|67[) and (|64[) we have 



V Ar+1 (R(X, Y;A 1 ,...,A r )Z) = A r+ ip~ x R{X, Y;A X ,...,A T )Z 

- RiAr+^X, Y;A±,..., A r )Z - R(X, Ar+ip^Y; A u ...,A r )Z 

r 

- R(X, Y;A U ..., ArXAr+ip^Z) - ^ R(X, Y;Ai,..., Ar+ip^A, A r )Z. 



(68) 



Combining the above we arrive at a version of (|66[) : 
R{X,Y;A 1 ,...,A r+1 )Z = 

= (l - |) {Ar+w^RtX, Y; Ai, . . . , A r )Z - R^^p^X, Y\ A\, . . . , A r )Z 

- R(X, A r+1 p- l Y; Ax,..., A r )Z - R(X, Y;Ai,..., A^Ar+ip' 1 Z) 

- Y; A\, . . . , Ar+ip^Ai, A r )Z^ 

- | {R{X, Y;Ai,..., Ar)Zp- x A r+r - R{Xp~ x A r +i, Y;A\,..., A r )Z 

- R(X, Yp^Ar+x; Ai, A r )Z - R(X, Y;Ai,..., A^Zp' 1 A r+1 ) 

- R(X, Y; A\, . . . , AiP^Ar+t, ...,A r )z\, 
i=i J 
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which is equivalent to 
R(X,Y;A 1 ,...,A r+1 )Z = 

= (1 - k) [Ar+xp-iRiX, Y; A\, . . . , A r )Z — R^+xP^X, Y ; Ax, . . . , A r )Z 

- R{X, Ar+xp-W; Ax,..., A r )Z - R(X, Y; A u . . . , A^Ar+ip' 1 Z) 

- J2 R(X, Y; Ai, A r+1 p~ 1 A i , A r )Z 1 
i=i J 

+ ^{p [p' 1 Ar+x^RiX, Y;Ax,..., A r )Z~\ 

- Riplp^Ar+x^^X^Y; Ax, . . .,A r )Z 

- RiX^lp^Ar+x^^Y]; Ax, A r )Z - R(X, Y;Ax,..., Ar)^ 1 Ar+x^ 1 Z}) 

- J2 R(X, Y;Ax,.. . .pIp-^+i.p- 1 ^], . . . , A r )Z 1 . 
»=i J 

(69) 

Now we can reverse the claim and using the exactly the same argument starting 
with the third equality in (|65|) we can prove that R(X, Y; Ax, . . . , A r )Z is the linear 
combination of terms Sp, where S is some word which is a product of the terms 
Xp^ 1 , Yp^ 1 , Axp^ 1 , . . . , Ap 1 r of the first order. Similarly we end up with 

R(X,Y;Ax,...,A r+1 )Z = 

= (1-k) {R(X,Y;Ai,...,A r )Zp- 1 A r+1 -R(Xp- l Ar+i,Y;Ai,...,A r )Z 

- R(X, Yp^Ar+x; Ax,..., A r )Z - R(X, Y; A u . . . , A r ){Zp- x A r+1 ) 

- J2 R{X, Y; Ax, A ip - l A r+ x, A r )Z 1 
»=i J 

+ | { \R{X, Y;Ax,..., A r )Zp- 1 ,A r+1 p- 1 ] p 

- RdXp^^r+xp-^p, Y; Ax, A r )Z 

- R(X, [Yp^^r+xp-^p; Ax,..., A r )Z - R(X, Y; Ax, . . . , A r )([Zp-\ Ar+xp^p) 

- J2 R(X Y;Ax,..., [A i p-\A r+ jp- l ]p, A r )Z 1 . 
»=i J 

(70) 

Now subtracting (|70|l from (|69[) and using the fact that 



[Ap 1 ,Bp 1 ]p = P [p 1 A,p 1 B]=-p\p 1 B,p 1 A] 
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for any A, B G H(n, C) and linearity of R(X, Y;At,..., A r )Z, we get that 
R{X, Y; At, A r+1 )Z - R(X, Y;At,..., A r+1 )Z = = 

= p [p^Ar+^p^RiX, Y;At,..., A r )zj - R(p[p~ 1 A r+ i,p~ 1 X],Y; At,..., A r )Z 

- R{X,p\ p - 1 A r+ i,p' 1 Y}-A ll A r )Z - R(X, Y; Ai, A r )(p[p _1 Ar + t,P^ Z]) 

r 

- R(X.,Y; Ai, . . . .pIp-^+i.p- 1 ^], . . .,A r )Z. 

■i=i 

(71) 

So in particular (|69[) is just 
R(X,Y;Ai,...,Ar +1 )Z = 

= (1 - k) j^+ip- 1 i?(X, K; Ai, . . . , A r )Z - RiAr+xp^X, Y;A U ..., A r )Z 

- R(X, Ar+tp^Y; At,. . .,A r )Z - R(X,Y; At, . . . , A r )(A r+1 p- x Z) 

- J2 R(X, Y;Ai,..., A r+ xP~ X Ai, . . . , A r )Z 1 
i=i ) 

= {1-k) P [p- 1 Ar+^RiX, Y; At, A r )Z - p^R^+tp^X, Y;At, A r )Z 

- p- x R{X, Ar+tp^Y; At,..., A r )Z - p- x R[X, Y; At, A^iA^p' 1 Z) 

-^p- x R{X, Y; At,..., A r+ ip~ 1 A i , A r )Z 1 . 
i=i ) 

Considering again (|64p and the first claim we get that the above is equivalent to 

R(X,Y; At,..., A r+1 )Z = (1 - k)D(R(X, Y; At, . . . , A r ))\p\ [A r+1 ], 

which is just (|63ll . □ 

Now again the Lie algebra f)p(V) is generated by the endomosphisms V r R. This 
means that the generated algebra grows as r increases and after some finitely many 
steps it stabilizes and taking higher covariant derivatives of R is unnecessary. Since 
the manifold P(n, C) is simply connected the holonomy group and the restricted 
holonomy group coincide, so f)p(V) = f) p (V). By the second formula in (|65|) and 
(|63|) we have the following 

Corollary 3 The Lie algebra f)p(V) is faithfully represented over the vector space 
V = H(n, C) (or V = H(n, R) ) with p : fj p (V) ^ End(V) given as 

p(W)Z=ZW + W*Z (72) 

for W G f)p(V) and Z £ H(n,C) (or Z G H(n,R)). 

We are in position to do a case by case analysis for different values of k. so(n, R) 
denotes the Lie algebra of skew-symmetric n-by-n matrices over the real field R, 
su(n, C) denotes the Lie algebra of skew-hermitian matrices with vanishing trace 
over C, s[(n, F) denotes the Lie algebra of traceless matrices over the field F. 
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Theorem 10 Let the smooth connected manifold P{n, C) with tangent space H(n, C) 
be equipped with real analytic connection 

V X Y P = DY[p][X p ] ~ | (Xpp^Yp + Ypp^Xp) (73) 

with k £ R. Then the holonomy algebra f)p(V) is as follows: 



P (V) 



the trivial algebra if k = 0, 2, 

su(n,C) «/« = !, (74) 

s[(n, C) e/se. 



/n i/ie case of the submanifold P(n, R) with tangent space H(n, R) we have 

Mv) 



t/ie trivial algebra if k = 0, 2, 

So(n,R) ifK=l, (75) 

s[(n,R) e/se. 



Proof By the conjugate invariancy of "H P (V) it is enough to consider the case when 
p = 7. 

Suppose k = 0, 2. Then the curvature ([62} of the connection vanishes, so f)p(V) 
is the trivial algebra. 

Suppose k = 1. Then the curvature (|62[) is nonzero, but is covariantly constant, 
all first and higher order covariant derivatives vanish due to Theorem [9] Therefore 
the manifold is a symmetric space that is very well known and the algebra f)p(V) 
by (|65p is generated by elements of the form [X, Y] where X, Y G H(n, F) . We have 
for all [X,Y] = W € f) P (V) that 

W* = [X,Y]* = -[X*,Y*] = -[X,Y] 

where * can be replaced by the transpose T over F = R. Also since TrW = 
Tr[X, Y] = we have f) p (V) = so(n, R) if F = R and f) p (V) = su(n, C) if F = C. 

Suppose k ^ 0,1, 2. Then by Theorem[9]the higher order covariant derivatives 
V 1 R (as we will see immediately) no longer vanish. Let W £ f)p(V). Then by (|65p . 
(j63|) and Corollary [3] 

^ = (f " T) (1 - «YD(b- l X,p- x Y\)\p)[A x , . . .,A r ], 

where X,Y € H(n, F). I.e. W is given by the linear combination of commutators of 
some n-by-n matrices over the field F, so TrW = 0. This tells us that 



f) P (V) C S [(n,F). (76) 

Now we will show that the generated algebra already stabilizes for r = 1. With- 
out loss of generality we can assume that p = I. Then we have to consider the 
generators of the form 



G = D^X^Y])^]^] = -[A X X,Y] - [X,A X Y]. 



(77) 
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Let 

E+ = \E ik + E ki if i^k, 
lk [Eii ifi = k, 

E ik =E ik - Eki 

where E^ is the matrix with zero entries excluding the (ik) entry which is 1. 
Then E% form a basis of H(n, R) and E~ k form a basis of the vector space of 
skew-hermitian matrices SH(n, R) over the real field R. The vector space SH(n, C) 
is defined similarly over C. Note also that H(n,C) H(n,R) SH(n, R) and that 
E tk E lm = in general. Suppose that A\ = E+, X = E+ y and Y = E+. Then by 

dZZI) v 

- -Eiz^ky^k + E ik E iz E ky - E ky h lz h lk + E lz b lk b ky . 
Using that E^ k = E^. and imposing restrictions z ^ k and y i we get that 

G= [ Ez v Hz^y, 
[Ezz-Eh if y = z. 

So the matrices G of this form span the whole s[(n,R), i.e. considering (|76p we 
have p)p(V) = sl(n, R) if F = R. Similarly if A x = X = £7^ and F = then 
we get the same generator G, so f)p(V) = sl(n, R) ©s[(n,R) due to E^~ k E^ n = if 
F = C; that is f)p(V) = s((n, C) in the complex case. □ 

By the proof of the previous Theorem 1 101 we see that Rj k i ;m = everywhere if 
and only if k = 0, 1, 2. This proves the following 

Corollary 4 The only matrix means which are affine means corresponding to sym- 
metric spaces are the arithmetic, harmonic and geometric means. 

Since we know the holonomy groups we can decide their metrizability. 

Corollary 5 The affine connections (|73|) are metric in the following cases: 

1. n = 1, 2, k arbitrary, F = R or C, 

2. n > 3, k = 0, 1, 2, F = R or C. 

Proof The case n = 1 is trivial. In [27] all irreducible holonomies of affine connec- 
tions are classified and metrizability is also dicussed. The metric connections were 
classified by Berger long ago. The holonomy sl(2, R) is isomorphic to sp(2, R) which 
is metric, there exists an invariant symplectic form. Also s((2, C) is isomorphic to 
so(l,3) which is also metric, there exists an invariant pseudo-Riemannian metric 
with signature (1,3). This isomorphic correspondence fails in higher dimensions 
n > 3, where the holonomies sl(n, F) (F = R or C) with representation over H(n, F) 
is non-metric. □ 

Remark 4 In the second case in Corollary [5] although there exists no metric struc- 
ture, however by inspecting the holonomy group T-L P (V) we get that there exist 
totally geodesic flat submanifolds. That is if we consider the subset D(n, F) of 
diagonal matrices of H(n, F) in both cases F = R or C, we get a totally geodesic 
Euclidean submanifold and a Riemannian metric on D(n, F) is given in the form 

Tr {p-^log^)} 

where log p (5) is the logarithm map given in (|37[). 
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So there exist no previously unknown affine matrix mean which correpsond to 
a Riemannian manifold. Although we have found a previously unknown, generally 
non-metrizable, one parameter family of affinely connected manifolds where the 
points of the geodesies are matrix means, in particular matrix power means. 
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